Are Gamma-Ray Bursts Standard Candles? 



o 
o 

(N 



O 

m 

Oh- 

6 



O 

o 
o 



x 



Li-Xin Li 

Max- Planck- Institut fur Astrophysik, 85741 Garching, Germany 



Abstract 

By dividing a sample of 48 long-duration gamma-ray bursts (GRBs) into four groups 
with redshift from low to high and fitting each group with the Amati relation 
log-Eiso = a + frlog-Epeak; I find that parameters a and b vary with the mean redshift 
of the GRBs in each group systematically and significantly. The results suggest that 
GRBs evolve strongly with the cosmic redshift and hence are not standard candles. 
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A remarkable achievement in the observation of GRBs has been the identifica- 



tion of several good correlations among GRB observables [see ISchaeferl (120071 ) 
for a review]. Based on several of those correlations, some people have eagerly 
proposed that GRBs are standard candles and can b e used to probe cosmology 
to very high redshift through the Hubble diagram (ISchaeferl. 120031 ; iDai et al. 



2004 : iGhirlanda et all 12004 : lLamb et all l2005t iFirmani et all 120061 ) . 



Type la supernovae (SNe la) have been considered as standard candles and 
applied to cosmology, which has produced a persuasive evidence that the Uni- 
verse is currently expanding with an accelerating speed (jClocchiatti 120061 and 
references therein). 

However, the situation of GRBs is very different from that of SNe la. The 
majority of the observed SNe la have redshift z ^ 0.1, so the evolution of SNe 
la plays a minor role. While for GRBs, all the correlations have been obtained 
by fitting a hybrid GRB sample with redshift spanning a very large range: 
from z ~ 0.1 up to z ~ 6. For objects in such a large range of redshift, it is 
hard to believe that evolution is not important. 

In addition, the physics of SNe la is much better understood than that of 
GRBs. We know that SNe la are produced by the thermal nuclear explosion 
of white dwarfs. For GRBs we are not aware of their progenitors, although 
many people think that long- duration GR Bs arise from the core-collapse of 
rapidly rotating massive stars (jPiranl . 12004 . and references therein). 
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The fact that people use a hybrid GRB sample to do statistics without dis- 
criminating the redshift distribution is of course caused by the fact that we 
do not have an enough number of GRBs with measured redshifts limited in a 
small range. Then, inevitably, the effect of the GRB evolution with redshift, 
and the selection effects, have been ignored. This raises an important question 
about whether the relations t hat people have f ound reflect the true physics of 
GRBs or they are superficial (IBand k, Preecd . 120051 ). 



A necessary condition for a class of objects to be standard candles is that they 
do not evolve with the cosmic redshift. Or, they evolve with the redshift but 
we know how they evolve (V. Petrosian, this proceeding). For GRBs, neither 
of these conditions is satisfied. 



I use the Amati relation as an example to test the cosmic evolution of GRBs. 
The Amati relation is a correlation between the isotropic-equivalent energy of 
l ong GRBs and the peak energy of their integrated spectra in the GRB frame 
flAmati et all . l2002f ) 
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With 41 long G RBs with firmly determined redshifts and peak spectral energy, 
Amatil (j2006al ) has obtained that a = —3.35 and b = 1.75 with the least squares 



method (E pea ^ in keV and E iso in 10 erg). 

If GRBs do not evolve with the cosmic redshift, we would expect that the 
Amati relation in equation (1) does not change with the redshift. Hence, a 
test on the variation of the Amati relation with the redshift could provide 
some constraint on the cosmic evolution of GRBs. 



For thi s purpose, I separate a sample of 48 long G RBs [41 GRBs f rom I Amati 
(J2006aJ), and seven additional Swift GRBs from I Amatil (j2006bl )] into four 
groups by the redshift, each group containing 12 GRBs: Group A — 0.1 < 
z < 0.84, (z) = 0.56; Group B - 0.84 < z < 1.3, (z) = 1.02; Group C - 
1.3 < z < 2.3, (z) = 1.76; Group D - 2.3 < z < 5.6, (z) = 3.40; where (z) is 
the mean redshift. 



Then, I fit each group by equation (1) and calculate the mean redshift, and 
check if the values of a and b evolve with the redshift. The results are shown 
in Fig. 1 (left panel; VL m = 0.3, fi A = 0.7, and H = 70 km s" 1 Mpc" 1 ). 

A least squares fit to the 48 GRBs as a single sample with equation (1) leads to 
a = —3.42 and b = 1.78, with the reduc ed \l = 5.9. The results are consistent 



with that obtained with 41 GRBs by lAmatil (2006a ). In the r ight panel of 
Fig. 1, I plot the deviation of fit [s; see iLi fc Paczyhskil (120061 )] against the 
mean redshift of GRBs. There is not a clear trend for s to vary with (z). But 
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Fig. 1. Left panel: Least squares fit to each of the four groups of GRBs by equation 
(1) (solid line). The two dashed lines mark the 1-a deviation of the fit. Right panel: 
The deviation of fit. Each circle corresponds to a group of GRBs. The star (s = 0.16) 
is the result obtained by fitting the whole sample (48 GRBs). 




01234 01234 



<z> (z> 

Fig. 2. The fitted value of a (left panel) and b (right panel) against the mean redshift 
of GRBs. Each data point with error bars represents a group of GRBs. The solid 
line is a least squares linear fit to a-(z) and b-(z). 

it appears that the deviation of fit of each group is smaller than that of the 
whole sample. This fact indicates that treating GRBs at different redshifts as 
a single sample may increase the data dispersion. 



I find that the values of a and b vary with the mean redshift of the GRBs 
monotonically. In Fig. 2, I plot a and b against (z). Clearly, a and b are 
correlated/anti-correlated with (z). The Pearson linear correlation coefficient 
between a and (z) is r(a, (z)) = 0.975, corresponding to a probability P = 
0.025 for a zero correlation. The correlation coefficient between b and (z) is 
r(b, (z)) = —0.960, corresponding to P — 0.040 for a zero correlation. 
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A least squares linear fit to a—(z) leads to 



a = -4.58 + 0.43 z , (2) 
with — 0.13. A least squares linear fit to b-(z) leads to 

b= 2.32 - 0.207 z , (3) 
with = 0.31. 

The results indicate that a and b strongly evolve with the cosmic redshift. 

I have used Monte-Carlo simulations to test if the variation of a and b with the 
redshift is caused by the selection effect. The results show that there is only a 
~ 4 percent chance that the observed variation is cau sed by the selection effect 



arising from the limit in the GRB fluence (ILil . 120071 ). Hence, the variation of 
the Amati relation with redshift that I have discovered may reflect the cosmic 
evolution of GRBs and indicates that GRBs are not standard candles. 

The results need to be tested with more GRBs. No matter what the conclusion 
will be (whether the variation of parameters is caused by the GRB evolution 
or by the selection effect), the results suggest that it is very risky to use GRBs 
with redshift s spanning a large range as a single sample to analyze physical 
relations among observables. Although I have only tested the Amati relation, it 
would be surprising if any of the other relations does not change with redshift. 
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